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Abstract 

A review is first presented of the Hall-Post inequalities relating TV-body to (TV — I)-body energies of quantum bound states. These inequalities 
are then applied to delimit, in the space of coupling constants, the domain of Borromean binding where a composite system is bound while 
smaller subsystems are unbound. 



L INTRODUCTION 

There are many examples, at various scales, of compos- 
ite systems at the edge between binding and non-binding. 
In nuclear physics, a proton-proton or neutron-neutron pair 
misses binding by a small margin, while a proton and a neu- 
tron form a rather weakly bound deuteron. The existence of a 
near-threshold state can induce dramatic consequences, for in- 
stance on fusion probabilities [ 1]. A pair of charmed mesons 
is presumably near the border separating stability from spon- 
taneous dissociation |2|. Atoms such as 4 Hc were for a long 
time believed to be unable to merge into a molecule. Recent 
studies indicates a tiny binding of the order of 1 mK for 4 He2- 
However, if one replaces one of the 4 He by an atom contain- 
ing the lighter isotope 3 He, then the 3 He 4 He is unbound. For 
a recent review on 3 HeAr 4 Hej\/ systems, see, e.g., Refs. QHl- 

An intriguing question is whether it is easier to bind three or 
more components than to form a mere two-body bound state. 
An answer is provided by the study of halo nuclei, which con- 
tain peripheral neutrons. Consider for instance the 6 He nu- 
cleus. It is stable against any dissociation, while the lighter 
5 He spontaneously decays into a neutron and a 4 He. In the 
(reasonable) approximation where the structure of the core is 
neglected, this means that the (a, n, n) three-body system is 
bound, while neither (a, n) nor (n, n) have a discrete spec- 
trum. 

This property of 3-body binding without 2-body binding 
was astutely named Borromean [5], after the Borromean rings, 
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FIG. 1 : Borromean rings 

'Dedicated to my colleague and friend Vladimir Belyaev at the occasion of 
his 70 th birthday 



which are interlaced in a subtle topological way (see Fig.^ 
such that if any one of them is removed, the two other become 
unlocked. The adjective Borromean is nowadays broadly ac- 
cepted in the field of quantum few-body systems. 

Borromean binding is intimately related to two other fas- 
cinating properties of few-body quantum systems. The Efi- 
mov effect |6] indicates that when the two-body energy van- 
ishes (e.g., by tuning the strength of the potential), a myriad 
of weakly-bound states show up in the three-body spectrum. 
This implies that the three-body ground-state already exists 
at this point. Slightly above the onset of two-body binding, 
the ratio E2/E3 of two-body to three-body binding energies 
is very small. By rescaling, one can reach a situation with a fi- 
nite 2-body energy, and a 3-body energy that becomes infinite 
when the range of the potential is made shorter and shorter: 
this is the Thomas collapse |7]. 

This review is organised as follows. In Sec.[n| the Hall-Post 
inequalities are briefly recalled. They are applied in Sec.lIIIIto 
constraint the domain of coupling constants leading to Bor- 
romean binding for bosons interacting through short-range 
forces. The difficulties arising in the case of fermions are de- 
scribed in Sec. II VI Borromean binding with Coulomb forces 
is the subject of Sec.|V] before the conclusions. 



II. HALL-POST INEQUALITIES 

A number of inequalities can be written down for binding 
energies in quantum mechanics if one splits the Hamiltonian 
into pieces (each piece being hermitian). Thus, for example, 

H = A + B + ■ ■ • => E(H) > E(A) + E(B) + ■■■ , (1) 

in an obvious notation where E(H) is the ground-state energy 
of H. Saturation is obtained if A, B, etc., reach their mini- 
mum simultaneously. If, for instance, H = p 2 — 1/r + r 2 /2 
describes the motion of a particle feeling both a Coulomb and 
an harmonic potential, then E(H) > (-1/2) + (3/2), cor- 
responding to an equal share of the kinetic energy. A slight 
improvement is obtained by writing H = [ap 2 — 1/r] + 
[(1 — a)p 2 + r 2 /2] , and optimising a. 

The reasoning can be applied to obtain a lower bound on 
3-body energies in terms of 2-body energies. This has been 
discovered independently by several authors working on the 
stability of matter 1 8] or baryon spectroscopy in simple quark 
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models 191. Let 



N 



(2) 



be the Hamiltonian describing a system of TV identical par- 
ticles interacting with pairwise forces, and E^(m,g) be its 
ground state energy . From the identity 



H N (m,g)=Y,H i ;l 1 [(N-l)m, 
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N -2 

i 

(i) 

where H N1 does not include particle i, one derives 
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N — 

> N{N ~ 1) E 2 (m(N-l),g)) 



(3) 



E N (m,v) > NE N -i (N - l)m 



(4) 



which can supplement any upper bound provided by a varia- 
tional method, to frame the exact energy. The generalisation 
to unequal masses or different potentials among the pairs is 
straightforward. For instance, 



H 3({mi};{gij}) 



pi 



Ami 



gi2V(r 12 ) 



(5) 



Numerical investigations show that the lower bound © is 
not very accurate and never reaches saturation. The relative 
wave function of particles 1 and 2 is of course modified by 
the presence of the third particle, and thus the three-body 
wave function does not describe optimally the (1,2) subsys- 
tem. However, this is a small effect which even disappears for 
the harmonic oscillator. The main source of inaccuracy in (0} 
is that the contribution of the (1,2) subsystem is replaced by 
the rest energy of an isolated (1,2) pair, although the (1,2) pair 
is not at rest in the whole system. 

To overcome this difficulty, a better decomposition was 
written down by Hall and Post 111 (il . and rediscovered in 
Ref. |11]. It involves the translation-invariant part Hn 



Hn — Hm — 



(Pi + Pi 



2Nm 



of each Hamiltonian Hm and reads 



H N (m,g)=J2 H 



N-l 



N(N - 2) g 



N-l ' N-2 



(6) 



(7) 



This leads to the new inequality 



E N (m,g) > E N -i m 



N(N — 2) 



> N ^ N -^E 2 



N - 1 
mN 



N 



(8) 



which is necessarily better than (0) since, mN/2 < m(N — 1) 
for N > 3, and for a given potential, the binding energy is a 



decreasing function of the constituent mass. This inequality 
is optimal in the sense that it can be saturated: this occurs for 
the harmonic oscillator and only in this case. The inequality 
also holds for the variational approximation E™ to the 
binding energy En, provided the sets of trial functions ip for 
N = 2 and N = 3 are consistent. In particular, one gets 
saturation, E™(m, v) = 3i?2 ar (3m/2, v) if each variational 
energy is calculated with a single harmonic-oscillator function 
ijj _/y oc exp(— ajv X) r 1j)- This is a rather frequent pattern 
that (with minimal restrictions) variational solutions share the 
properties of the exact ones. For the virial theorem this was 
underlined by Fock and rediscovered by many authors 1 12]. 

The case of particles with unequal masses has revealed 
some surprises. The simple extension (here for N — 3) of 
the decomposition 



H 3 [m,i;Vi 



13 i ■ 



Oi\2 



mip 2 — m 2 Pi 
mi + m 2 



r, 
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(9) 

does not always give saturation for the harmonic oscillator. 
The remedy was found in Ref. iTPjll for N = 3 particles and 
extended in Ref. for N — 4: one should introduce more 
freedom in the decomposition of the Hamiltonian and find the 
optimal choice by suitable adjustment. 

For N = 3, the decomposition is written as 



H 3 = BP + a 12 



Pi ~ x 3 P 2 

l + X 3 



+ gVi 



12 



(10) 



where P — J^Pi- The identification gives the irrelevant vec- 
tor B and the inverse masses in terms of the free param- 
eters Xi. The lower bound on E 3 is expressed as a sum of 
2-body energies. It depends on these X{ and is maximised by 
varying these parameters, for each given set of potentials Vij . 
In the case of the harmonic oscillator with unequal masses 
and even unequal strength factors, saturation is reached, i.e., 
the exact energy always coincides with the lower bound. 

For iV = 4, one should go a step further, and write the 
decomposition as 



H A = B.P + a 12 



,X\2,iPi 



9V1 
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(11) 



The momentum ^\ xi 2 , , p t is constrained by translation in- 
variance, and 2:12.2 — ^12.1 = 1, to make it conjugate to 
r 2 — n, This leaves two free parameters for each pair, which 
are adjusted to optimise the lower bound. 

Refs. fill [l3l [A contain several examples, where the 
three- or four-body energy is computed and compared to the 
above lower bounds. For confining potentials, one is not too 
far from the limiting case of the harmonic oscillator for which 
the inequality is saturated, and the inequalities ^ provide 
a very good approximation to the exact energy. The quality 
of the approximation deteriorates for short-range or Coulomb 
forces. 
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III. BORROMEAN BINDING OF BOSONS 



IV. THE CASE OF FERMIONS 



The Hall-Post decomposition was used in the previous sec- 
tion to obtain inequalities on energies at given coupling con- 
stant g. They can also be used to provide relations between 
coupling constants at given energy E. In particular for E — 0, 
one can derive inequalities on the critical couplings g^, where 
gN is the minimal strength required to achieve TV-body bind- 
ing. 

In one- or two-dimensional quantum mechanics, any attrac- 
tive potential supports at least one bound state |15]. More 
precisely, g2 = if the integral J d n xV(x) is negative. For 
n = 3 dimensions, the problem is more subtle 1 16]. An attrac- 
tive potential with short range, say gV(r), requires a minimal 
strength g > g 2 to achieve binding of two constituents of unit 
mass, and mg > g 2 for constituents with mass m. A clas- 
sic paper by Blatt and Jackson 1 17] gives estimates of g 2 for 
simple potentials such as Yukawa. 

In this context, the phenomenon of Borromean binding is 
expressed by the property that the critical coupling 33 required 
to bind three bosons is smaller than g 2 . More generally, g^ < 
gN-i- Several questions can be raised: what is the typical 
range of values for (73/(72, 94/93, etc.? Are there rigorous 
constraints on these ratios? Are these constraints saturated for 
particular potentials? What can be said in the limit of large 
TV? 

Some results are given in Refs. flsLflslEoTl . which contain 
references to earlier papers. In particular, the decomposition 
Hn <x Hn-i implies that Hn hardly explores the do- 
main of binding with negative expectation values, if all Hn-i 
remain positive. Thus 

Ng N > (TV ~l)g N -i > ...2g 2 . (12) 

For simple monotonic potentials, one gets 33/32 — 0.80, 
well above the rigorous bound 33/52 > 2/3. Still, there is 
a remarkable window of about 20% in the coupling constant, 
where three-body systems are bound without two-body bind- 
ing. One also gets 34/32 — 0.67, typically, which, when com- 
pared to 33/32 — 0.80, reveals a window of about 13% for 
four-body binding without three-body binding. 

For bosons, the situation is rather simple. If a potential does 
not succeed in binding a pair of bosons, but is predominantly 
attractive, that is to say has a negative scattering length, it will 
bind a sufficiently large number of bosons. 

We have seen that Hall-Post inequalities are saturated for 
harmonic potentials. A pure V(r) oc r 2 does not give Bor- 
romean binding, as gV would confine even for 3 — > 0+. But 
a potential like V oc r 2 — C in the domain of the relevant wave 
functions, and vanishing at very large distances r, would give 
33/32 ~ 2/3, i.e., nearly saturate the bound [20]. This corre- 
sponds to an attractive potential with an external barrier. 

More important for physical applications are potentials with 
an internal repulsive core. Their wave functions significantly 
depart from those of the harmonic oscillator and the ratio 
33/32 become closer to 1. Examples are given in 1 20] . 



In general, an inequality such as (|8}, though remaining 
valid, is of little interest for fermions. The ground state E 2 
corresponds to a symmetric orbital wave function, with, for 
instance, the spins arranged in a singlet state to fulfil the an- 
tisymmetry requirement. On the other hand, three or more 
electrons can never be in a fully-symmetric orbital wave func- 
tion. Thus the actual En, TV > 3, will be an excited level of 
the Em entering the inequality (|8j. 

A minimal knowledge of the structure of the wave function 
is thus needed to write down an useful inequality. Consider 
for instance three spin 1/2 fermions in a symmetric spin (or 
spin-isospin, or spin-isospin-colour) state. Then the space 
wave function of the ground state should be a L p = 1+ an- 
tisymmetric wave function, a prototype being the harmonic 
oscillator state p x Aexp[— a(p 2 + A 2 )] in terms of Ja- 
cobi coordinates p and A. The two-body subsystems are 
in an antisymmetric state with orbital momentum £ = 1. 
Then Ez(m,g) > 3E 2 (3m/2, 3) results from the identity 
i?3(m, 3) = ^2 H 2 (3m/2, 3) applied in the subspace of anti- 
symmetric wave functions. 

For more than TV = 3 fermions, or other spin states of TV = 
3, the problem is more complicated but still, some rigorous 
results can be obtained. 

The basis is the decomposition of TV-particle representa- 
tions of the permutation group in terms of TV — 1 states [21]. 
For instance, in a spin 5=1/2 state of three electrons a given 
pair is half of the time in a singlet state, and half of time in a 
triplet state. This means a proper bound on E3 involves, for 
E 2 , and average over the £ = and the I = 1 ground state. 
For other cases, one needs some group theoretical calculations 
l2lll . This gives nice local constraints, i.e., links from En to 
En -1 (see also the improvements proposed in Ref. l22lo . But 
this approach fails when extrapolated from E 2 to the large-TV 
limit of En- This means that much progress remains to be 
made. Checking the constraints with exactly solvable poten- 
tials 1 23J gives valuable insight into the problem. 

Note that the question of writing inequalities for fermion 
energies is rather old. In the course of the pioneering papers 
on the stability of matter, a very clever decomposition was 
written down 

hn+i =y^jy 1 

i 

a a<(3 

where the superscript in X^ 1 ' indicates that particle i is omit- 

(i) 

ted. Each h J is an independent-particle Hamiltonian with 
ground-state energy yV €j, in terms of one-particle energies 
ej, properly piled up as explained in textbooks of elementary 
chemistry (2 particles in Is, 6 in 2p, etc.), to ensure antisym- 
metry. However, the recoil of, e.g., the {2,3 . . . , TV} subsys- 
tem against the first particle is not taken into account. As a 
consequence, the inequality is never saturated. 

Significant progress was achieved by Basdevant and Martin 
12411 . who used sophisticated convexity inequalities to derive a 
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bound for power-law potentials, that becomes saturated in the 
limit of harmonic forces. 



V. BORROMEAN MOLECULES 

Borromean binding as described in Sec. [Ill] seems at first 
hardly conceivable for a Coulomb potential, since changing 
the strength by an overall factor g simply results in a mere 
rescaling by g 2 of all binding energies E^. However, if the 
stability of ions and molecules does not depend on the overall 
strength of the interaction, it relies on an adequate balance be- 
tween attraction and repulsion, and, for a given set of charges, 
on the ratios of the constituent masses which are involved. 

Consider for instance a set of three masses {rrii} carry- 
ing charges {qi} = ±{+1,-1,-1}, (times an overall fac- 
tor). Systems with identical or nearly identical inverse masses 
m^ 1 and m^ 1 , such as H^ (e~,p,p), Ps - (e + ,e~,e~) or 
H~ (p, e _ ,e _ ), or neighbouring configurations, are stable, 
while less symmetric systems such as (e~ , p, e + ) or (p, p, e~) 
spontaneously decay into an atom and an isolated charge 1 25] . 

Note that the stable systems with three unit charges are not 
Borromean, since there are always two stable atoms among 
the subsystems, {1,2} and {1, 3} in our notation. This means 
that one can construct each stable system step by step: two 
charges form a stable atom, which attracts the third charge. 

For N = 4 unit charges {qi} = +1, +1, —1, —1, there are 
new possibilities. The most familiar cases are the following: 

• unstable configurations, such as (p. e + . p. e~) l3ll . 
which spontaneously splits into a protonium and a 
positronium, 

• stable systems, the simplest being the Ps2 molecule 
(e + , e + , e~, e~), whose all three-body subsystems 
±(e + ,e~,e _ ) are stable, 

• less stable systems, such as the positronium hydride 
(p, e + , e~ , e~ ) , which can be viewed as a stable H~ fix- 
ing a positron, or a stable Ps~ attached to a proton, or, 
more interestingly, as an unstable (e~, p, e + ) stabilised 
by the addition of a second electron. 

More recently, a new category of four unit-charge systems was 
revealed 12611 . They are stable, but all three- body subsystems 
are unstable. Hence, they are Borromean, in the sense that 



they cannot be built by adding the constituents one by one. An 
example is (p, d, p, d), with a proton, a deuteron and their an- 
tiparticles. Indeed, a study by Mitroy Il27il . who uses Varga's 
stochastic variational method |28], indicates that asymmetric 
ions (m + , M~ , m~) are stable in the range 

M 

0-70 < — < 1.69 , (14) 

TO 

while other studies, by Bressanini et al. 1 29], who use a dif- 
fusion Monte-Carlo approach, or Varga 1 30 ] with a stochastic 
search of the range coefficients of a Gaussian parametrisation, 
have established that (M + , to + , M~, to - ) remain stable for 



Then, for M/m ~ 2, and for neighbouring configurations 
(M + , to + , M'~, to' - ), the four-body system is stable while 
none of its three-body subsystems is stable. 



VI. OUTLOOK 

In molecular, nuclear, or hadronic physics, there are fragile 
and subtle quantum states, whose stability relies on the joint 
efforts of many subsystems, which by themselves do not con- 
tain enough attraction to be stable. It is of course tempting 
to extend to political or sociological situations: l'union fait la 
force, as it is said in a popular French proverb. 
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